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Abstract 

We applied finite difference time domain (FDTD) algorithm to the study of field and intensity 
correlations in random media. Close to the onset of Anderson localization, we observe deviations 
of the correlation functions, in both shape and magnitude, from those predicted by the diffusion 
theory. Physical implications of the observed phenomena are discussed. 

PACS numbers: 
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Light propagating in a random medium undergoes multiple scattering. As a result, it 
forms a complex interference (speckle) pattern when it emerges from the medium, which can 

nn 

be described statistically|l|,|2j. There exist longer range correlations, that can be detected as 



correlations between two distant speckles 



^. These longer range correlations in transmit- 



ted intensities received much attention in context of mesoscopic electronic systems jsl, where 
they result in universal conductance fluctuations. For electronic systems only limited infor- 
mation (fluctuation of conductance) can usually be deduced from experiments. For hght, 
however, momentum- (angular) and spatial-resolved measurements are possible. The latter 
yield much more detailed information about the system. The connection between electrons 
and photons is established by identifying the Landauer conductance g with the total trans- 
mission coefficient (summed over all incoming and outgoing channels) ^. The majority of 



experimental and theoretical efforts (see Refs. 



8| and references therein) to date 



are concentrated on systems in the regime of diffusive transport, g ^ 1. The intensity 



correlation function (ICF) is defined as 



C(Ar Au) = inr + ^r,u + Au)Iir,u)) _ 
^ ' ^ (/(r + Ar,z/ + Az/))(/(r,z/)) ' 



(1) 



where r and u are spatial coordinate and frequency respectively. When length L of t 
random medium is greater that its width W, C is invariant with respect t o r, and isotropic 



for Ar, as long as one avoids the evanescent zone on the output surface ly] . Theoretically, 



he 



9 



based on pairing of incoming and outgoing channels, three contributions to ICF have been 



identified 
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8[: short-range Ci, long-range C2 and infinite-range C3. Deep into diffusion 



regime (yf ^ 1 in a wave-guide geometry 
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8| Ci ^ \CE{Ar,Au)\^ ~ 1, C2 ~ l/g and 



C3 ~ 1/(7^, making the values of C2 and C3 small. Here, 



CE{Ar,Au) 



{E{r + Ar, + Au)E*{t, v)) 
(/(r + Ar, z/ + Az/)) V2(/(r, z/))i/2 ' 



(2) 



is the field correlation function (FCF). Essentially, the nonpertubative nature of the crossover 
from diffusion to localization limits the amount of information that can be obtained 



theoretically 



ll| . The purpose of this letter is to introduce a numerical method for studying 



correlation functions (CFs) based on FDTD 



12l | algorithm. This method allows us to solve 



Maxwell's equations for electromagnetic field at every spatial grid at every time step. It 
makes no assumption about the scattering strength, accounts for all interference phenom- 
ena, and makes it possible to obtain CFs in both diffusion and localization regimes. Here 
we apply our method to systems with values of g from 4.5 to 1. 

We consider a 2D system as shown in inset of Fig. 1: parallel plate metallic waveguide 
filled with circular dielectric scatterers of refractive index n = 2 and diameter d = lA cm. 



We choose our parameters close to the microwave experiments of Ref. 1^. The scatterers 
were randomly positioned (without overlapping) with a fixed filling fraction / = 0.3. Metallic 
walls of the waveguide insured 99.9% reflectivity, and uniaxial perfectly matched layer 12 1 
boundary conditions were applied outside the waveguide. A TM polarized broadband pulse 
with center wavelength A = 2 cm (close to one of Mie resonances of the scatterers) was 
launched via a point source placed one wavelength away from the left surface of the random 
medium inside the waveguide. A temporal discrete Fourier transform (DFT) was applied to 
electric fields observed at a series of vertical planes separated by A/2 at the output end of the 
waveguide. By virtue of DFT, we obtained the continuous wave response of the system for 
a number of wavelengths Aj within a 4% range of A. This range was significantly narrower 
than (i) the width of the broadband pulse, and (ii) the range in which physical parameters 
of the system start to deviate from those at A. The distance between two consecutive Aj was 
chosen to be equal to the average mode spacing. Spatial field and intensity CFs are obtained 
by setting Au = in Eqs. I1I2I We averaged over 48 configurations and 60 Aj, similar to 



the microwave experiments of Ref. [13|. In addition, due to the invariance of correlations in 
the empty part of the waveguide, we averaged over 10 observation planes. Up to 100 spatial 
points equally spaced near the center of each plane were sampled depending on Ar. For 
spectral CFs, Eqs. (Ill2j) with Ar = 0, the averaging procedure is similar except that the 
number of Aj sampled depends on the value of Az/. 

Waveguide geometry makes our system quasi-lD, where localization length C, oc Nl. N 
is the number of waveguide channels, and / is the transport mean free path. With an 
increase of L, the system crosses over from diffusion to localization while / is kept constant. 
According to the diffusion theory the dimensionless conductance of a quasi- ID system g = 



a qi : 



{-K /2)n~^jjNl / L\ where L' = L + 2zb accounts for the boundary effect jl4|. the extrapolation 



length Zh is usually of order Z, rigjj = c/ve-, ve is the energy transport velocity In a 
passive system g is equal to the Thouless number which is defined as the ratio of the diffusion 



mode linewidth to average mode spacing 
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In our system, Z = 1.8 cm, rS^jj = 1.77(see 



below). The width of the waveguide, W = 20 cm {N = 2W/X = 20) and lengths L = 20, 40 
and 80 cm yielded the values of g = 4.4,2.2 and 1.1 respectively, and allowed us to study 
the correlation functions near the onset of localization. 

Zh can be determined by fitting the real part of spatial FCF with the formula 
^ ,A ^ vr(V0'^o(A;Ar) + 2sin(ZcAr)//i;Ar 

Ce{M = /, , ^ > (3) 

where k = 2tt/X, and Jo is Bessel function of zero order. The imaginary part of C£;(Ar) 



should vanish due to isotropy 



9|, which is confirmed by our calculation in which its value was 



ess than 10 . We derived the above expression in 2D following the 3D derivation of Ref. 



17l |. Eq. (jni) gives an excellent fit for all systems studied with Zb/l = 0.8 in Fig. 1. Absence 
of the deviation of C£;(Ar) from the expression derived in diffusion regime becomes apparent 
once one recognizes that Ce{^t) contains only information about short range correlations. 



It reflects the correlations on length scale of /, which is much smaller than C,- n^eff 
expression for g can be calculated n^jj = (1 x Wair + n x Wgcat) / {Wair + Wscat)-, where 
WaiT and Wscat are the energy stored in air and scatterers respectively. Wair and Wscat were 
determined numerically to give riflj ~ 1.77. The physical parameter yet to be obtained is 
the transport mean free path /. This can be done by fitting the spectral CF to the complex 



function 



18 
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C^(Az/) = gL'/sinhgL', (4) 



with q = y TiAu/ D{l—i). I enters Eq. Q through L' and the diffusion coefficient D = fE//2. 
Fig. 2 shows the fitting of the real and imaginary parts of C£;(Az/), from which we find 
I = 1.8 cm. The half width at half maximum (HWHM) of |C£;(Az/)p should coincide with 
1.46 times the diffusion mode linewidth Su = D / L' , at least in the diffusion regime. The 
inset in Fig. 2 shows that the obtained / and are consistent. A slight deviation in the 
system of g = 4.4 is attributed to its short length. 

The spatial dependence of long-range contribution to ICF was derived in diffusion up to 



1/g'^ order [13 



C(Ar)-|C.(Ar)|-=l^ + ^l (5) 



4 8 \ 1 + \CE{Ar) 



Fig. 3 shows both magnitude and normalized profile of Eq. (0). The numerically calculated 
C(Ar) — |C£;(Ar)|^, shown in Fig. 3, reveals that its spatial profile is independent of 
g. Specifically, the ratio of its value at Ar = to that at Ar oo remains equal to 
two. Stronger oscillations in the numerical data are likely due to the finite width of the 
waveguide. However, the magnitude of the long-range contribution is significantly enhanced 
due to localization effects, far beyond the diffusion prediction up to order l/g"^. 

In Fig. 4, as g decreases, the long range contribution to spectral ICF, C(Az/) — |Cs(Az/)|^, 
is broadened when Az/ is normalized to Su. We ascribe this effect to strong fluctuations close 



to the localization threshold. Namely, a few (more conducting) channels with larger than 
average linewidth dominate ICF, leading to its spectral broadening. 

In conclusion, using the FDTD algorithm, we studied field and intensity correlation func- 
tions close to the onset of localization. In this regime neither experiments nor analytical 
theories have given such detailed information about correlation of intensities transmitted 
through a random medium. 
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Figure captions 

Figure 1. Real part of spatial FCF. Squares, circles and triangles correspond to the systems 
with g = 4.4, 2.2, and 1.1 respectively. Solid curve represents the fit of Eq. Q with 
Zb/l = 0.8. The inset shows the geometry of our system. 

Figure 2. Real (empty symbols) and imaginary (full symbols) parts of frequency FCF. 
Solid and dashed curves represent real and imaginary parts of C^; given by Eq. Q 
with / = 1.8 cm. The inset compares found from HWHM of |C£;(Az/)p to D/L'^. 
Both quantities are normalized to average mode spacing. Symbol notations are the 
same as in Fig. 1. 

Figure 3. The magnitude of long-range contribution to ICF versus dimensionless conduc- 
tance g. Symbol notations are the same as in Fig. 1. Solid line represents diffusion 
expansion formula Eq. (jSJ at Ar = 0. The inset shows long-range contribution to 
ICF normalized to one at Ar = 0. Solid, dashed and dot-dashed curves correspond to 
samples with g = 4.4, 2.2, and 1.1 respectively. Thick solid curve plots Eq. (0). 

Figure 4. Frequency dependence of long-range contribution to ICF normalized to value at 
Au = . Symbol notations are the same as in Fig. 1. 
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